An optimized three-dimensional time-space domain staggered-grid finite-difference method

Numerical simulation of three-dimensional (3D) seismic wavefields forms the basis research on migration methods 3D data based wave equations. Because precision wavefield extrapolation determines imaging accuracy to a certain extent, it is very important study how enhance forward modeling wavefields. Thus, we build an optimized staggered-grid finite-difference (SFD) method with high two-dimensional (2D) modeling. Since generates corresponding difference coefficients by utilizing least square (LS) minimize objective function constructed time-space domain dispersion relation acoustic equation, our LS-based SFD can effectively in theory compared Taylor-series expansion (TE), especially for large wavenumber range. Examining numerical dispersion, algorithm stability and computational cost, compare three conventional TE-based illustrate demonstrate its effectiveness feasibility. The examples from different models suggest that generate less higher than other methods, but condition stricter cost slightly higher.